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Abstract. We prove a formula expressing the gradient of the 
phase function of a function / : R'^ H> C as a normahzed first fre- 
quency moment of the Wigner distribution for fixed time. The for- 
— mula holds when / is the Fourier transform of a distribution of com- 

r^ , pact support, or when / belongs to a Sobolev space i?''/^+^+^(]R'') 

where e > 0. The restriction of the Wigner distribution to fixed 
^ , time is well defined provided a certain condition on its wave front 

set is satisfied. Therefore we first study the wave front set of the 
V^ , Wigner distribution of a tempered distribution. 

< 

^ . 0. Introduction 

-(— > 

This paper treats a time-frequency version of the following trivial 

observation in Fourier analysis. Let f{t) = Ce^'^*^°'*, C G C \ 0, ^o ^ 

M'^, be a nonzero complex multiple of a character on M'^. Its Fourier 

'pi . transform is / = 06^^ so the frequency ^o may be expressed using the 

'nI" ! Fourier transform as the normalized first order moment formula 

l> 

q: (0.) .„^|| 

o 

O : where (/, is the vector (/, ^ = ((/, Q)'^^^ G R'^ and ^j : M'^ ^ M is 

coordinate function j, I < j < d. 

We will deduce a time-frequency version of this formula for more 
^ . general functions, which looks like 

^Varg/(t) = |^Ml|) WGM'^: /(t)^0, f e ^^'{R'). 
In the formula (10.21) Wf denotes the Wigner distribution, defined by 

for / G ^(M'^). For functions in ^S"{M.'^) which are not multiples 
of characters e^'^*^"'*, the frequency is not well defined. Therefore it is 
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replaced in fl0.2p by the natural generalization 

^Varg/(t) = i-(9,arg/(t))^t,, 

that is, a normalized gradient of the phase function. We use the term 
instantaneous frequency, taken from the engineering literature [Hll] , for 
this quantity. Thus (10. 2p may be seen as a time- frequency version of the 
observation flO.ll) . We shall also prove a version of (10.21) for functions 
/ that belong to a Sobolev space iJ'^/^+^+^(M'^) where e > 0. Then the 
distribution actions {Wf(t,-),^) and (VF/(t, ■),!) may be replaced by 
Lebesgue integrals. 

In order to prove (10. 2p we need to restrict the Wigner distribution as 
Wf f-^ Wf{t, ■) to fixed time t e R'^. For / e y"(R'^), the restriction is a 
map ^'(M^^) H-)- ^'(M'^). Restriction of a distribution to a submanifold 
is possible under certain conditions on the wave front set p]. More 
precisely, the restriction defines a well defined distribution provided 
the normal bundle of the submanifold has empty intersection with the 
wave front set of the distribution. Thus we are led to study the wave 
front set of the Wigner distribution first. We pursue this study in 
somewhat greater generality than actually needed in order to prove 
formula ([021) • 

We define the space WFW-^ of tempered distributions such that the 
wave front set of the Wigner distribution is directed purely in the fre- 
quency direction, and the space WFW^ of tempered distributions such 
that the wave front set of the Wigner distribution is nowhere parallel to 
the time direction. The latter space admits restriction Wf i-)- Wf(t, ■). 
We show the inclusions C^^ ^ WFW^ and C,^„ C V^^n ^ WFW^. 
Here C°^^^ denotes the space of smooth functions on M.'^ such that a 
derivative of any order is bounded by a constant times a fixed polyno- 
mial. C^^ contains ^S". The space Kon is the subspace of ^'{R'^) 
such that the short-time Fourier transform (STFT) V^f satisfies 

sup {{x,r]))"'\V^f{x,r])\ < oo Vn > 0, 

\r)\>B\x\ 

for some B > 0. This means that the STFT decays polynomially in a 
conic neighborhood of the frequency axis. 

Recently Guo, Molahajloo and Wong have studied the instantaneous 
frequency and its relation to the modified Stockwell transform [7J. For 
other aspects of wave front sets and time-frequency analysis we refer 
to the recent papers [2[l3tfTH[T2] . 

1. Preliminaries 

The Schwartz space ^(M'^) consists of smooth functions such that 
a derivative of any order multiplied by any polynomial is uniformly 
bounded. Its topological dual y'(M.'^) is the space of tempered dis- 
tributions. We denote by C^(]R'^) the space of smooth and compactly 
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supported functions, and ^'{R'^) is its topological dual, the space of dis- 
tributions. The compactly supported distributions are denoted S"{R'^). 
The normalization of the Fourier transform for functions / G =5^(]R'^) 
used in this paper is 



^/(o = m = / /(x)e-2--«d 






where x ■ ^ denotes the inner product on M.'^. The inverse Fourier 
transform is then 



The Fourier transform extends by duality to a homeomorphism on 
^'(M.'^). For s G M the Sobolev space H'^iM.'^) is defined as the subspace 
of / G y(R^) such that / G Ll^iR"^) and 

1/2 



H^ = (^Jjo'vm'd^^ <oo 



where (0 = (l + |eP)'/'. 

Translation is denoted by (Tyf){x) = f{x — y) and modulation 

by M^f{x) = e^'^*^'^/(x) for functions of one M'^ variable, and by 
{Ty,^f){x,z) = f{x-y,z-w), M^J{x,y) = e2-(-«+^-')/(x,i/), re- 
spectively, for functions of two M'^ variables. The short-time Fourier 
transform (STFT) of / G J^'{M.'^) with respect to a window function 
^ G y(R'^) [51E] is defined by 

where (■, ■) denotes the conjugate (for consistency with the L^-product) 
linear action of ^'{R'^) on C^iR"^), or y{R'^) acting on S^{R'^). Thus 
the inner product on L^, also denoted by (•,■), is conjugate linear in the 
second variable. We denote the linear (without conjugation) action of 
distributions on test functions by (■, ■), and thus {u,ip) = {u,Ip). 

We use the symbol C for a positive constant that may change value 
over inequalities and equalities. The space Co(M'^) consists of con- 
tinuous functions that vanish at infinity. Thus / G Co(M'^) means 
that for any e > there exists a compact set K^ C R'^ such that 
X ^ K^ ^ \f{x)\ < e. The symbol C^iW'-) is the space of functions 
such that all partial derivatives of order not greater than k are contin- 
uous everywhere, and C°°(M'^) = {^^^^C^{R'^). 

We recall the definition of the C°° wave front set of m G ^'(M*^) [5l|8], 
denoted WF{u). Let F C M*^ \ denote an open conic subset, where 
conic means ^ G F ^ a^ G F for all a > 0. The wave front set is 



4 P. BOGGIATTO, A. OLIARO, AND P. WAHLBERG 

defined as the complement 
WF{u) = {R'^ X M<^ \ 0) \ 

{(x,0: 3V^GCr(M'^):^(x)^0,3rcM^\0: ^ e T, 
snp{ri)"'\%jju{ri)\ < oo Vn > 

We have PiWF{u) = singsupp(M) where Pi denotes the projection on 
the first M"^ variable. The singular support singsupp(M) is the comple- 
ment of the largest open set where u is C°° . 
The cross-Wigner distribution is defined by 



,, ,, Wat,0= /(t + r/2Mt - r/2)e~2--W 

(1.1) JRd 

where 

K,{x,y) = {x + y/2,x-y/2) 

and ^2 denotes partial Fourier transformation in the second M"^ vari- 
able. The Wigner distribution of a single function is denoted by Wf = 
Wfj. The original definition (11. ip for f,g ^ ^{R'^) extends to f,g G 
y\R'^), and then Wf^g e y'{R'^'^). 

The Wigner distribution has been studied thoroughly from many 
points of view, in particular quantum mechanics ^SE\, pseudo-differential 
calculus |S] and signal analysis [6j. In signal analysis the Wigner dis- 
tribution has been studied as a way to represent signals simultaneously 
in the time and frequency variables. The Wigner distribution satisfies 
many properties requested by an ideal time-frequency representation, 
among which the most important include the following, which holds 
for example when / G S^{R'^). 

(1.2) WM.T.f = T,,r,Wf, x,r]e M^ 



:i.3) // Wf{t,Odtd^=\\ 

:i-4) / Wf{t,i)di = \f{t)\ 
:i-5) / Wf{t,odt = \m\ 
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(1.6) ^Varg/(t) = ^r^Jf/V?5 Vt G M'^ : /(t) ^ 0. 

Formula (II. 2p says that the Wigner distribution respects time- frequency 
shifts, (II. 3p says that its integral equals the squared energy of the func- 
tion, and (11.41) . (11.51) say that Wf has the correct marginal properties. 
This admits the interpretation of Wf as a distribution of the energy 
of / over (t, ^) G M^'^. Alternatively, in quantum mechanics, Wf is 
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interpreted as a phase space probability density. However, these in- 
terpretations are in general not possible, since Wf is not nonnegative 
everywhere unless / is a generalized Gaussian of the form 

(1.7) f{t)=exp{-TTt- At + 2TTb-t + c) 

where c e C, b e C^, A e C'^'^ is invertible and Re A > 0. This is 
Hudson's theorem [5l|6]. The function / has Wigner distribution 

Wf{t,^) = Cexp (-27r (t • Re At - 2Re bt)) 

X exp (-27r(^ + Im At - Im b) ■ (Re A)"^(^ + Im At - Im b))) 

and 

— V arg /(t) = -Im At + Im b. 

27r 

Thus Wf is a function mainly concentrated along the submanifold 

{(t,(27r)-Varg/(t)) : teW'jCR''''. 

The formula (II. 6p expresses a generalization of this observation, in 
the sense that the instantaneous frequency is a normalized first order 
frequency moment of the Wigner distribution, for fixed t G M"^. It 
is well known in the applied literature |I] where it is derived without 
precise assumptions ioi d = 1. 

Janssen ^^0\ has studied the question whether the Wigner distri- 
bution (for d = 1) may be concentrated on a curve in the phase space. 
Under some assumptions on the curve it turns out that it must be a 
straight line and the distribution / is either a multiple of a Dirac distri- 
bution or a degenerate Gaussian of the form (11.71) with Re A = 0. This 
means that Wf is of the form Wf{t, = C6o{^ - {27i)-^V arg /(t)) (as- 
suming 6 = 0), i.e. supported on the subspace {(t, (27r)~^Varg/(t)) : 
t e M} C M^. For other functions the Wigner distribution gives a dis- 
persion. However, (II. 6p shows that the mean value of Wf{t, ■) for fixed 
t e M"^ agrees with intuition since it equals the instantaneous frequency. 

This paper concerns assumptions that imply that (II. 6p holds true. It 
is relatively straightforward to relax / G .^{W^) into / G iJ'^/^+^+^(]R'^) 
where e > (see Proposition 14. ip . since Wf is then a continuous func- 
tion which admits restriction Wf i— >■ Wf{t, ■) without problem. As 
another assumption we use / G ^S"{R'^) (see Proposition 14. 2p . which 
implies that Wf G ^'(M^'^). This case is more subtle since restriction of 
a distribution to a submanifold is not always possible. However, there 
are conditions on the wave front set that are sufficient for a restric- 
tion to be well defined and continuous. Hence we need to investigate 
the wave front set of the Wigner distribution. This problem is a gen- 
eralization of the study of the singular support of Wf, considered by 
Janssen ^^. 

The paper is organized as follows. We investigate the wave front 
set of the Wigner distribution in Section |2j In Section [3] we study 
sufficient conditions for the restriction operator Wf H- Wf{t, ■) to be 
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well defined for t G M'' fixed. Finally in Section H] we prove formula 
dnSD for / G ^^'(M'^) and for / G i/^/^+i+e^^d) ^j^g^g ^ > g. 

2. The wave front set of the Wigner distribution 

Let / G o$^'(]R"'). Two natural questions are to compare singsupp(/) 
and Pi singsupp(PV/), and to compare WF{f) and Pi^2,WF{Wf). Here 
Pi : M^'' I—)- M'^ denotes the projection on the first variable, Pi{t, ^) = t, 
t, ^ G M*^, and Pi^s : M^'^ i— )■ M.^^ denotes the projection on the first and 
third variables, Pi^3{t,^;ri,x) = {t,ri), t,^,ri,x G M.'^. 

Let / = 1. Then singsupp(/) = and WF{f) = 0. Moreover, 
Wf = 10 5o so singsupp(Vr/) = M'^ x and WF{Wf) = {{t,0;0,x) : 
teR'^, X G M"' \ 0} (see Example Q. Hence Pi singsupp(iy/) = M'^, 
Pi,3WF(Wf) = M'^ X 0. Thus 

Pisingsupp(iy/) ^ singsupp(/), 

Pi,3WF{Wf) ^ WF{f), f G y'{R''). 

We do not know whether any of the following inclusions hold for 

/ G y{R'^). 

(2.1) singsupp(/) C P^ singsupp(iy/), 

(2.2) WF{f) C Pi^^WF{Wf). 

Note that the inclusion (12. 2 p is stronger than (12.11) . In fact, assume 
([22D. We have sing supp(iy/) = Pi,2WF{Wf) and therefore PiWF{Wf) 
Pi singsupp(iy/). The assumption implies 

singsupp(/) = PiWF{f) C PiPi^3WF{Wf) = PiWF{Wf) 

= Pi sing supp (Wf), 

and thus ([2l]) follows from (12:21) . 

To produce a counterexample to the inclusion (12. ip . it suffices to find 
a function / which is not C°° at t = whose Wigner distribution Wf 
is C°° in a neighborhood of (0,^) for all ^ G M'^. There is some weak 
evidence that such a function may exist: in fact the Wigner distribution 
is regularizing in the sense that 

(2.3) W : L^R"^) x L\R'^) ^ CoiR^"^) 

continuously. To see this, let (/„) be a sequence of functions /„ G 
C^{R'^) such that ||/-/„||i2 ^ as n ^ oo. Since Wf„ G ^(M^^) and 
since ||W/,c,||l°° < 2*^11 /||i2|| 51 1|2,2 by the Cauchy-Schwarz inequality, 
we see that Wf is the uniform limit of Co(M^'^) functions, and since 
Co(M^'^) equipped with the supremum norm is a Banach space we have 
Wf G Co(K^'^)- 

A characterization of the wave front set WF{f) in terms of the as- 
ymptotic behavior of Wf for large frequencies is given in Corol- 
lary 3.28]. Let (f> G ^(M'^) be even, nonzero and define the dilation 
0^(a;) = s'^/^cpis^/'^x) for s > 0. Then (^0,^0) ^ WF{f) if and only if 
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there exists a neighborhood U of (to,^o); conic in the second variable, 
such that for any a > 1 

(2.4) sup s''\{Wf*W^s){t,sO\<oo Vn>l. 

{t,0&U, a''^<\^\<a, s>l 

This criterion says that the convolution Wf * W^s decreases faster 
than any polynomial in the frequency direction in a neighborhood of 
(to, ^o); conic in the second variable. Note that the function W^s (t, ^) = 
W(j,{s^^'^t,s~^^'^^), with which Wf is convolved, concentrates around 
zero in the t variable and spreads out in the ^ variable as s — t- +oo. 

In particular we have that / is smooth in a neighborhood of to if and 
only if there exists a neighborhood V of to such that for any a > 1 

sup s''\{Wf*W^s){t,s^)\<oo Vn>l, 

teV, a-i<|g|<a, s>l 

which means that Wf * W^s decreases faster than any polynomial in 
any frequency direction in a neighborhood of to. 

The criterion (12. 4 p says that the wave front set WF{f) can be char- 
acterized by Wf. Roughly speaking, the microregularity of / is char- 
acterized by the asymptotic behavior of Wf at infinity in conic fre- 
quency domains. (More precisely, this behavior concerns the convolu- 
tion Wf * W^s and not Wf.) 

In the remainder of the this section we will investigate the wave front 
set of the Wigner distribution Wf for / G ^'(M*^). We introduce the 
following two subspaces of ^'(M'^) for this purpose. Here we under- 
stand by subspace a subset which is not necessarily lineaiu. 

Definition 2.1. 

WPW^iR"^) = {f e y{R'^) : WF{Wf) C M^^ x (0 x M^ \ 0)}. 

Thus WFW^{M.'^) consists of the tempered distributions such that 
the wave front set of the Wigner distribution is directed purely in the 
frequency direction (or is empty). 

Definition 2.2. 

WFW^iR"^) = {f e y{R'^) : WF{Wf) n M^^ x (M'^ \ x 0) = 0}. 

The definition says that WFW^(R'^) consists of the tempered dis- 
tributions such that the wave front set of the Wigner distribution does 
not contain vectors purely in the time direction. Obviously 

WFW^iR'^) C WFW^iR'^). 

The following example shows that for /(t) = exp(7r2t ■ At), A sym- 
metric, we have 

f eWFW^{R'^)\WFW^(R'^), AeR'^'"^\0. 



We do not know whether WFW^ or WFW^ are hnear subspaces of y". 
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Example 2.3. Let f{t) = exp{7Tit ■ At) where A G M''^'^ is a symmetric 
matrix (/ is sometimes called a chirp [B]). Then 

Wf{t,0= I exp(-27r2r-(e-At))rfr = 5o(e-At), 

singsupp(Vr/) = {{t,At), t e M"'}. 

Since the transformation T : M^'^ i-> M'^, T(t, ^) = ,^ — At has surjective 
differential, the distribution Wf = 6q o T e ^'(M^"') is well defined 
[H Theorem 6.1.2]. Let ip E C^(R^'^) and let (ipn) C C^(R'^) be a 
sequence converging in ^'(M'^) to (5o. The continuity statement of [H 
Theorem 6.1.2] gives 

^Wf{r], x) = {So o T, M_n-^(f) = lim (V'n o T, M_^ _^(/)) 

n— >-oo 

= lim // ^„(e-At)(^(t,Oe-'"'(''-*+"-«)rftde 



r( 



where x(^) = V'l^, ^t) G C° 

Let (t, ^) G singsupp(Vl//), i.e. ^ = At and suppose ip{t,At) ^ 0. If 
?7o + Axo 7^ then there is a conic neighborhood F containing (r/oj^o) 
such that \ri + Ax| > C > when |(77,x)| = 1 and (r7,x) G F. This 
gives \ri + y4x| > C|(?7,x)| for {r],x) G F and thus, using the fact that 
XG^(M'^), 

sup \{r],x)\^\ipWf{r],x)\<Cn sup [(r/, x)|"|?7 + Ax]"" < oo 
(T?,x)er ('7,3;)er 

for any n > 0, so (t, At; ?7o, Xq) ^ WF{Wf). 

On the other hand, we may use the following result [S] for u G 2i'{^^) 
and v?i e C;?°(]R'^). If ^j{x) ^ for all j > 1, supp((/9j) -^ {x}, and 
^jM does not decay polynomially in any conical neighborhood of ^ for 
any j, then (x, ^) G WF{u). 

Let v? satisfy V9(t, At) ^ and y? > 0. Then £(0) = / </)(«, Au)du > 

and (y5pVj(— Ax, x) = £(0) 7^ for any x G M^. Therefore ipWf is not 
polynomially decreasing in any conic neighborhood of the manifold 
{{i],x) : rj + Ax = 0}. By shrinking the support of ip we may thus 
conclude 

WF{Wf) = {(t. At; -Ax, x), teR, X G M \ 0}. 

For A = we thus have / G WFW^{R'^) but if A ^ then / G 
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We are interested in WFW^{R'^) since / G WFW^{R'^) is a suffi- 
cient condition for the restriction operator Wf H- Wf(t, ■) to be well 
defined from ^'(M^^) to ^'{R'^) for all t e R'^ (see Section ED- We 
are interested in WFW^iW^) since it is easier to prove inclusions of 
familiar spaces in WFW-^{R'^). 

Next we define a linear space of smooth functions, whose derivatives 
are slowly increasing, uniformly with respect to the order of derivation. 

Definition 2.4. The space Cg^^(]R'^) consists of smooth functions such 
that a derivative of any order a e N*^ is bounded by a constant C^ times 
a fixed polynomial, that is / G Cg^^(M'^) if there exists A^ > such 
that 

(2.5) |<9"/(x)| < Ca{x)^ Vx G M^ C„ > 0, a G N'^. 

Note that C3^„(M'^) C Oij{R'^), where Of^{R'^) is the space of smooth 
functions such that the derivatives satisfy (12.51) where N may depend 
on a. For more information on (9^(]R'^) we refer to [13] . 



Lemma 2.5. 



^^'(R'') C C^JR''). 



Proof. Let / G ,^(g"(R'^) which means that there exists u G S"(R'^) 
such that f{x) = {u,(f)x) where ipxiO = exp(27rzx ■ ^) [H Theorem 
7.1.14]. Denote by A^ the finite order of the distribution u G (S"{R'^). 
By [HI Theorem 2.1.3] we have d°'f{x) = {u, d^cjyx) for any a G N'^. For 
some compact set K cR'^ containing the support of u, this yields 

\d"f{x)\ = \{u,d:<Px)\ <c^Y. sup|a,-af0.(oi 



<C^{xy\ aeW 



u 



Remark 2.6. If / = 1 = ^5^ then Wf = I ® 6o, and WF{Wf) = 
(R"^ X 0) X (0 X M'^ \ 0) according to Example O This shows that the 
wave front set of the Wigner distribution of a function / G ^(S"(R'^) 
in general is nonempty. 

Remark 2.7. Note that the inclusion ^S"(R'^) C C^^^(M'^) is strict. 
For instance C^(R'^) C C3^„(M'^) \ ^(^'(R'^). 

Next we shall make some preparations in order to introduce a linear 
space which is larger than C^^^(M'^). 

Lemma 2.8. If f e ^'(R'^) and ip G ^(M"^) then 



WfW^iv, x) = VJ{-x/2, T^l2)VJ{xl2, -^2), ?], x G 
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Proof. We compute 

= e"*''-(^+") {if® if ok) ( ^^, y-z + x 



e 



-'^<y+^)^{y + x/2)^{z-x/2) 



= M^/2T^x/2Viy)M^^/2T^/2(p{z). 
This gives, using [8l Lemma 7.4.1] and the fact that W^ is real- valued 

= {f®JoK,^i\M,,,W^)) 
= {f^7,^2'\M,,,W^)oK~') 



n 



= (/ (g) /, M^/2T_^/2V O M_r,/2T^/2(p) 

= VJi-x/2,r]/2)V^fix/2,-r]/2). 

Lemma 2.9. Let g E C°°(R.^'^) and suppose 

\git, 1 < C((t, 0)^ for some N > and all t,^ E M^ 
r = {(t,0 : 1^1 > C\t\} for some C > 0, and 

sup {{t,0)"]g{t,0\<oo \fn>0. 
(t,«)6r 

IfipE ^(M^rf) then we have for any C > C and T' = {(t,0 : l-^l > 

cm 

(2.6) sup {{t,0)^\9*^{t,0\<oo Vn>0. 

(t,?)er' 

Proof. We have 

1^ * ^{t, 1 < / / l^(^ -x,^-ri) 1 1 v?(x, ri)\dxdri 



+ \9{t - x,^ - ri)\\ip{x,r])\dxdr) 

Consider the first integral Ji. For any C > C and the corresponding 
cone r = {(t,0 : |^| > C'\t\} C T we have: (t,0 e T' and ((x,r/)) < 
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((t, CjY^"^ imply that (t— x, ^—77) G V provided that ((t, ^)) is sufficiently 
large. In fact, these assumptions imply 



and the quotient approaches one as ((t,0) ^ °^5 because |^| > C'\t\. 
Thus we have for any integer n > 



h<C^ {it-x,^-r]))-^\ipix,r])\dxdv 

(2.7) 



J jR2d 

provided ((t, 0) is sufficiently large. 

Next let us look at the second integral l2- We have for any L > 

I2<Cl If {{t-x,^-r])f{ix,r]))-'^dxdr] 

J J {{x,n))>{(t,0)^/^ 

<CL{{t,Of 11 {{x,v)f-'^dxdv 

= CL((t,0)^ [[ (1 + |a:p + |r^|2)(^-^)/2rfa;rfr/ 

(2.8) J J {{x,r!))>{iWy/^ 

= CL{{t,0)'' I (1 + r2)(^-^)/V^-irfr 

ir>(((i,0>-l)l/2 

provided that L > N + 2d and ((t, ,^)) is sufficiently large. Since L > 
is arbitrary, (12.71) and (12. 8p prove (12. 6p . that is g*{p(t,^) decays rapidly 
(polynomially) for (t,^) G T' . D 

Remark 2.10. Obviously Lemma 1^^ is invariant under a change of roles 
of the variables, that is, with cones F, F' of the form F = {(t, ^) : |t| > 

cm,c>o. 

We are now prepared to define the following linear subspace oi S^'{M.'^) 

Definition 2.11. Let y? G ^(M"^) \ 0. 
(2.9) 

{f eS^'iR'^): 3B>0: sup {{x,7])y\V^ f{x,r])\ < 00 Vn > 0}. 

\v\>B\x\ 

According to this definition V^^{M'^) consists of tempered distribu- 
tions such that the STFT decays rapidly in some conical neighborhood 
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of the frequency axis of the form {{x,t]) G M^'^ : |?7| > -B|x|}. Ob- 
viously V^j^iM.'^) is a hnear space. A priori V^^lW^) depends on the 
window function (f G ^(M'^), but the next lemma shows that this is in 
fact not the case. 

Lemma 2.12. If f e Kon(^"') then f G V^t^(R'^) for any ip G y(R'^) \ 
0. 

Proof. Let / G V^^{W^). According to [6l Lemma 11.3.3] we have 

\V^f{t,0\ < C\VJ\ * \V^^\{t,0, t,C e M'^. 

Furthermore, we have V^(p G ^(M^^)^ y^f G C°°{R'^'^) and by [6l 
Theorem 11.2.3] there exist C, A^ > such that 



The result now follows from Lemma 12.91 D 



As a consequence of the latter lemma we may denote V^^„(]R'^) = 
Vcon(I^°'), which is understood to be defined by an arbitrary ip G ^{R'^)\ 
0. Another consequence is that if (12. 9p holds, that is V^f{x,ri) decays 
rapidly in a conical neighborhood of the frequency axis {{x,ri) G M^*^ : 
|?7| > i?|a;|}, then rapid decay holds for a neighborhood of the form 
{{x,!]) G R^'^ : |?7| > B'\x - y\}, B' > B, for any fixed y G R'^. In fact 
we have VT_yipf{x, rj) = V^f{x — y, rj) which gives 



sup ((x,r7))"|1/^/(x,r/)| 

{x ,ri):\ri\> B' \x~y\ 

= sup {{x-y,r]) + {y,0)y''\VTy^f{x-y,ri)\ 

\v\>B'\x-y\ 

<C{yr sup {{x-y,r])r\VTy^f{x-y,r])\<oo Vn > 

\ri\>B'\x—y\ 



by the proof of Lemma 12.121 and Lemma 12.91 
Proposition 2.13. We have the inclusions 

(2.10) C's'£w(^'^) ^ WFW^iR'^), 

(2.11) CZJ^'') C Kon(M'^) C WFW^il 
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Proof. Let / G (^.^^(M'^) and cp G y{R'^). For a G N"^ integration by 
parts and (I2.5p give for some A^ > 



\v''V^f{x,r])\ 



77^ jR'i 



(2.12) <^\^n. f /A^/+-^\~^-d-^ 






/3<a 

Let 5 > be arbitrary and define the cone F = {{x,rj) G M.^'^ : |?7| > 
i?|x|}. Then if (x, rj) G F, we have by (12.121) for any n > 

which leads to 

sup ((x,r7))"|V;,/(x,r/)| < C sup {ri)''\V^f{x,ri)\ < oo 

for any n> Q. Thus Vipf{x,ri) decays rapidly in a cone \r]\ > B\x\ for 
any B > 0. This proves the first inclusion in (12. lip . 
By Lemma [2.81 we have 

(2.13) sup {{x,ri)y\W/W^{ri,x)\ <oo Vn > 

\v\>B\x\ 

for any 5 > 0. Let x e C^(M?'^) and ^ G ^(M'^) satisfy (xW^</.)(i,0 ^ 
for arbitrary fixed (t,0 e M^^. We have xW^<^ e (7^°^(M2'i), and 

(2.14) =F(iy/W^x)(^,a:) = W^^ * xiv,^)- 

We have W^^ = Wj *W^ e C°°(M2rf) and W/W^ has polynomial 
growth by [H Theorem 7.19] since 1?/ G ^'(M^"!) and W^^ G ^(M^'^). 
This means that there exist C,N > such that 

|Vr^^(r7,x)| <C((r/,x))^ r/, x G M''. 

Now fl213|) . (127141) . Lemma ES] and Remark [2:T0] says that for any 
B' > B and corresponding cone F' = {\ri\ > B'\x\} C F, it holds that 
^{WfWipx){Vy ^) decays rapidly in the cone F'. Hence ^{WfW^px){j], x) 
decays rapidly in a cone |?7| > B\x\ for any i? > 0. Therefore / G 
WFW^{W^). The inclusion fl^H?]) is therefore proved. 

It remains to prove the second inclusion in (12.111) . Let / G Kon(IR°')- 
This assumption and Lemma [2l8] imply that WfWip{ri, x) decays rapidly 
in a conic neighborhood |?7| > B\x\ for some B > for any (f G 
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y{R'^). Hence by Lemma Ell ^{WfW^x){v,x) = W^^ * x{v,x) 
decays rapidly in a conical neighborhood of the form \ri\ > B'\x\ where 
B' > 0, for any x e C^{R'^'^). This means that (t,^;r/,0) ^ WF{Wf) 
for any t, ^ G M'^ and r/ G M'^ \ 0, that is / G WFW^{R'^). D 

Remark 2.14. A combination of Lemma [23] and Proposition l2.13l fl2.10p . 
gives ^(f'(]R'^) C WFW-^(R'^). This admits the interpretation that the 
very high (analytic) regularity of / G ^S''{M.'^) is reflected in the fact 
that Wf inherits smoothness in the time direction. In fact the wave 
front set of Wf is directed purely in the frequency direction and has no 
component in the time direction. 

Remark 2.15. We note that an alternative proof of the inclusion 

can be deduced from the results by Toft [171[T8]. In fact, [HI Proposi- 
tion 1.8 and Theorem 4.1] imply the following result. 

If / G ^'(M*^) and (0, 0; r], x) i WF{Af) 

then (t, e; r/, 27rx) i WFi^^f) for all (t,0 G M^'^. 

(Note the factor l-n that appears in front of x, due to different normal- 
izations of the Wigner distribution.) Here Af is the so called ambiguity 
junction of / [5,^, normalized for / G S^{R^') as 

Af{^, X) = (-) f fit - r^)JitT^e'''-^dt. 

One can show with computations resembling those of the proof of 
Lemma 12.81 that 



AfA^{r],x) = \V^f{TTx,-TTr])\^, t],x e 



^d 



for / G y(R'^), (fi G y(R'^). Here (^(t) = y?(-t). Using Lemma 
12.91 as in the proof of Proposition 12.131 with (xA^)(0,0) 7^ we ob- 
tain {0,0;ri,0) i WF{Af) for r/ G M'^ \ provided / G Kon(M'^). An 
application of Toft's results finally gives Kon(M'^) ^ WFW^i^^). 

We may summarize the inclusions we have found as follows. 

^S'{^^) C C,7„„(M'^) C Kon(M'^) C WFW+i^^), 

C,^„(M"') C WFW^i^^) C WFW*{^^). 

The next example shows that Kon(K^) ^ WFW^i^'^) and (^.^^(M'^) C 

Example 2.16. Let r^o G M \ and set /(t) = exp(27r2r/ot2). Then 
by Example O we have WFiy^j) = {{t,2r]ot; -2r]ox,x), t eM.'^, x G 
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R'^ \ 0}, so / ^ WFW^{R'^). If we choose ip{x) = exp(-27rx2) then it 
can be verified that 

\Kf{x,v)\ =C'»?oexp (- 2 iVox-vf^f] , Cvo > 0- 

If i? > 0, \r]\ > B\x\ and < e < 1/2 we have 
V 



VOX-- 



> Y - mm > |r/| I - - — ) > e\r]\ 



provided B is sufficiently large. Thus 

sup {{x,r])r\VJ{x,r])\<C sup {vrKf{x,r])\ 

\v\>B\x\ \ri\>B\x\ 

27r£2 



< C sup {?])"' exp ( — - — — ^I^P ) < C)0, 

which means that / G Vcon(IR'')- Hence we have shown 14on(I^'^) ^ 
WFW^{R'^). For the function / we also have / ^ C,^w(^'^) which 
shows that C^^^lM"^) C KonlK"')- In fact, we have 



''slowv 



where Pa is a polynomial of order |a|. If we suppose that / G C^ '^'^^ 



then (12 .Sp gives for some A^ > 

sup |a"/(t)|(t)-^ = sup b„(t)|(i)-^ < oo 



for all a G N'^, which is a contradiction. Thus / ^ C^^„(]R'^). We note 
furthermore that a cone \ri\ > B\x\ where Vipf{x,ri) decays rapidly has 
B > 2|r7o|. Therefore it is not always the case that the cone of decay 
for elements in V^on(I^'^) can be arbitrarily large, that is, B arbitrarily 
small. 

Let A G W^^'^ be symmetric. The transformation Ta defined by 

TAfit) = exp{mt ■ At) fit) 

is continuous on ^(M'^) and extends to continuous transformation on 
e5^'(M'^). The Wigner distribution is transformed according to 

(2.15) WT,f{t,0 = Wfit,^-At). 

The next result treats the invariance and noninvariance under Ta of 
four of the spaces introduced above. However, for Kon(I^'^) "we cannot 
prove or disprove invariance. 

Proposition 2.17. Let A G R'^^'^ be symmetric. 

(i) ^(^'iR'^), C,^„(M°') and WTW^iR'^) are not invariant under Ta- 
(ii) WTW^(R'^) is invariant under Ta- 
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Proof, (i) Since 1 G ^(^'(R'^) C Q^^(M'^) C WTW^(R'^) by Lemma 
[231 and Proposition 1^1^ it suffices to show that g ^ WTW'^iR'^) 
where g{t) = exp('7rit ■ At). According to Example 12.31 we have 

WF{Wg) = {(t, At; -Ax, x) : teR'^, xeR'^\ 0}. 

Picking X eR'^ such that Ax ^ reveals that g ^ WTW-^{R'^). 

(ii) According to (1^13]) we have Wr^f = Wf o Q ioi f e y(R'^) 
where Q denotes the invertible matrix 

/ \ _ Tr^2dx2d 



«-' -A I 

According to 0, Theorem 8.2.4] we have 
WF{Wf o Q) 

= {{t, e; Q\v, x)) : (g(t, 0; V, x) e WFiWj), t, ^,v,xe R'} 
= {{t,^;r]- Ax, x) : {t,^- At; r], x) G WF{Wf), t, i,n,x e R''}. 

If / G WFW^{R'^) and {t,i;r],x) G WF{Wf) then x ^ 0, which 
implies that 

WFiWT^f) pl (M^"' X (M^ \ X 0)) = 0. 

Hence T^/ G WFW^{R'^). D 

2.1. An elaboration of the inclusion fl2.10p . In this subsection we 
prove a result related to the inclusion C3^„(M'^) C WFW^{R'^). We 
will prove a stronger statement under a stronger hypothesis. The con- 
clusion again includes the fact that the wave front set is directed purely 
in the frequency direction. On top of that we add the statement that 
the location of the singular support is included in M*^ x 0. 

We begin by recalling the definition of two classes of symbols: the 
H5rmander classes ^"^(XxM^) [8] and the Shubin classes S'^{R'^) [16]. 

Definition 2.18. Let X C M^ be open, m e R, p e (0,1], and 5 G 
[0, 1). Then ^'"^(A: x M^) is the subspace of a G C°°(X x M^) such that 
for every a G N'^, /3 G N" and every compact A^ C X there exists a 
constant Cm,a,(i,K > so that 

is satisfied for 2; G A', r G M'^. 

Definition 2.19. For m G M, < p < 1, T'^{R'^) is the subspace 

of / G C°°(]R'^) such that for every 7 G N'' there exists a constant 
Cm,7 > so that 

\dif{t)\<Cm,,{tr-'\'\ 

is satisfied for every t (iR"^. 



THE WAVE FRONT SET OF THE WIGNER DISTRIBUTION 17 

Remark 2.20. S^^iX x R'^) and r™(M'^) are Frechet spaces with re- 
spect to the best constants appearing in the estimates. If we extend 
Definition 12.191 to p = we note that 

r-(M'^) c cz^R'^) = U r-(M'^), < p < 1. 

We set r = M" x [Rf \ 0). A phase function on F is a real- valued 
smooth function that satisfies the conditions: 

(i) ^{z, \t) = X4>{z, t) for {z, r) e F, A > 0. 

(ii) V0(z, t) 7^ for every {z, r) G F. 
We recall the meaning of oscillatory integrals of the type 






where a G S'^g{R"'~^'^) and is a phase function on F (cf. [HI Theorem 

7.8.2]). 

Proposition 2.21. (i) For fixed u G C^{R^), and fixed phase function 
(j), the map defined by the absolutely convergent integral 

(2.16) C^(R"+^) 3 a— ¥ I^{u) = f e'^^''^^a{z, t)u{z) dz dr 

has a unique extension to a continuous functional on S^^{R^^'^) for 
?7i G M, p G (0, 1], (5 G [0,1). Hence I^{u) is well defined for every 
u G C^{R'') and a G ^;^5(M"+'^). 

(ii) For fixed a G ^"^(M""'"'^) the map (well-defined from (i)): 

(2.17) Cr(M") 3u-^ i;{u) 

is a distribution in ^'(M") of finite order, which is indicated by 

n(.)= I e"^^-'^'^a(;T)dT 



and called an oscillatory integral of symbol a and phase 0. 

A general result on the wave front set of oscillatory integrals is [8|, 
Theorem 8.1.9], which follows. 

Proposition 2.22. For the distribution /? we have the inclusion 

(2.18) WF{i;) C {{z, V.0(z, r)) : V^z, r) = 0}. 

Now we use the previous result in the study of the wave front set of 
the Wigner distribution. More precisely, in the following proposition 
we prove that, considering a function / in a Shubin class, the wave 
front set of the Wigner distribution Wf is not only "vertical" in the 
dual variables {ri,x), but the singular support of Wf is also contained 
in the "horizontal" subspace ^ = of the space variables (t,^). 
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Proposition 2.23. If p > and f e r™(M'^) then 

(2.19) WF{Wf) C {{t,0;0,x) G M^"' : t G M^a; G M"' \ 0}. 

Proof. If / G r;^(M'^) we have / ® / o k G ^^^(M^d). in fact, for every 
a, /3 G N'^ and for t in a compact set K C M'^, we have: 

\d?d^f{t + T/2)fit-r/2)\ 

= d^ Yl C/^,/^'^^'/(t + T/2)d^-^'fit - r/2) 

/3'</3 



E E Ca,a'C/5,/3'5"'+^7(t + T/2)d---'+^-^'f{t - r/2) 

a'<a/3'</3 

< C E E ^^ + r/2)™-^l"'+'^'l (t - r/2)"^-/'l"-"'+/^-/3'l 

a'<Q/3'</3 
^ (-;( Y^ V^ /^\|m-p|a'+/3'|| + |m-p|Q-Q'+/3-/3'|| / \m-p|Q'+/3'|+m-p|a-Q'+/3-/3'| 

a'<o/3'</3 

where Petree's inequahty {t ± r/2)'' < Cs(t)'''l(r)*, s G M, has been 
used. 

We set now z = {t, .^) G Mf x M^ and remark that a symbol a(t, r) G 
^^^(Mf X M^) can also be seen as a symbol in S'^^siK'^ x K), and 
likewise a phase function 0(t, r) on M^ x M^ is also a phase function on 
M.'jf' X M^. We can therefore apply Proposition 12.221 in the case n = 2d 
with 



ait,^,T) := /(t + r/2)/(t-r/2) G S.^l^lM?,^^ x R^) 
and phase function on MH x M.J 

0(t,e,r) = -27^^r. 
Thus we have Wf = I^, and 

V,0 = (0, -27rr) G M^^ 
V,0 = -2< G M^ 

which means that fl2.18p implies the inclusion fl2.19p . D 

3. Restriction of the Wigner distribution to fixed time 

In this section we shall study the restriction operator of a distribution 
F G ^'(M^o!) to the submanifold t x M'^ C M^d fo^ f ^ ^d g^^^^ ^j^j^^ 

is denoted 

(3.1) RtF = F{t,-). 
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This map is not well defined for any F G 3!'{E?'^). But according 
to [SI Corollary 8.2.7], the restriction (13.1 p gives a well defined element 
in ^'{W^) provided WF{F) n A^"* = 0, where 

(3.2) Nt = {tx M^) X (M*^ X 0) 

is the normal bundle of the submanifold t x M^ C M^'^. For / G 
WFW^{W^) we have WF{Wf)r]Nt = for any t e M"', so [8, Corollary 
8.2.7] implies that the restriction Wf{t,-) gives an element in ^'(M'^) 
for any t eW^. 

The distribution f = Sq = 6q{R'^) has Wigner distribution Wf = 
^0 ® 1 [S], which has wave front set 

WF{Wf) = (0 X R"^) X (M^ \ X 0). 

This example shows that WF(Wf) fl A't = is not always satisfied. 

We will study the continuity properties of the restriction (13. ip , and 
for that purpose we need the following definitions and results from [8]. 
For a closed set F C M'^ x (M'^ \ 0), conic in the second variable, we 
define 

^^{W^) = {ue ^'{R'^), WF{u) C F}. 

If u,Uj G ^f(M'^) for 1 < J < oo we say that Uj -)■ m in ^f(M'^) 
provided 

(3.3) Uj — ^ M in ^'(M'^), and 

(3.4) snp\^n^{^{u-uj)m\^0, J^oo, Vn G N, 

for a closed conic set V CR^ and ip G C^iW^) such that 

F n (supp (^ X y) = 0. 

It is sometimes more convenient to use instead of (13. 4p the equivalent 
requirement [8] 



1^^' sup sup |^ri^(OI < oo Vn G N. 

For a closed set F C M^'^ x (M^^ \ 0), conic in the second variable, we 
set 

Tt = {(e, x) G R^' : 3r/ G M"^ : (t, e; r/, a;) G F}. 
We will use Theorem 8.2.4]. Let F C M^d x (M2d\o) be a closed set, 
conic in the second variable, such that F fl A^^ = 0. Then [SI Theorem 
8.2.4] implies in particular that Rt is a continuous map 

(3.5) Rt:^r(^^'^)^^r^(^'^)- 

Next we study F G C^^(M^°') and the following two operators: 
(i) the restriction ([31]) to the submanifold t x M'^ C M^d f^^ ^ ^ j^d 
fixed, and (ii) Fourier transformation in the second variable F i— >■ ^2F. 
The following results say that these two operators commute for F G 
Cg^^(M^'^). First we need a lemma. 
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Lemma 3.1. Let F e CZ^iR'^'^), T = M^d x f/ where 

(3.6) f/ = {(r/,x) gM^^\0: |r/| < |x|}, 

let X e C^iR^), x>0, x{x) = 1 for \x\ < 1, Xj{x) = xi^/j), J > 
integer, and Fj = Fxj ® Xj ■ Then we have 

(3.7) ^2Fj — > ^2F m %{R?'^) as j -^ cx). 

Proof. If ^2V e C,°°(M2^) then ip G ^(M^'^) and we have, since Fj -^ F 
in ^'(M^rf) as J ^ oo 

i^2Fj, ^2^) = {Fj, if) ^ (F, if) = {^2F, .^2^), J ^ oo. 

Thus ^2Fj -^ ^2F in ^'{R'^'^) as j -^ oo and the first criterion fl3^ 
for (13.71) is proved. To prove the second criterion we use ([S3D'. Since 
'^'^gi-V - ■)) = M^v9 for 9 ^ '^iM'^) we have for ^ e C^{R'^'^) 
^{ip^2F,)ir],x) = ^-\ip^2Fj)i-v,-x) 

= ^i^Fj * ^~V(-^, -x) 

= (^fiF„^-V(-^- 

= (F,,^fi(^-V(-r7 

= (F„M_,,o^2-V(-,-a:--)) 

F{t, r)e-2-*-'^^2-V(t, -X - T)xj{t)Xj{r)dtdT 



-X--)) 



ff Fit, r)e-2-^*-^^2"V(t, -X - T)xAr)dtdT, 

J jR2d 



for i > J for sufficiently large J . In fact, the last equality holds for such 
j since ^2^Lp{-,y) has support in a fixed compact set, independent of 
y G R'^. The assumption implies that 

(3.8) |5fF(t,r)|<C,(t)^(r)^ a G N^ 

for some Cq,, A^ > 0. By means of integration by parts, (13. 8p and the 
observation that ^2^'{) G S^ , we obtain for j > J and any a G N'^ 

\v''^{^^2F,)iv,x)\ 

27r^)-l°l9r (e-^-^*-") F{t, T)^2~Mt, -x - r)xM)dtdT 

dtdr 



< 



(5<a 
I3<a 

<Co.Y.^x) 



Y^cJf d^Fit,T) dr'^2-Mt,-^-r 

(t)^(r)^(t)-^-'^-i(r + xy^'-'^-^dtdr 
{t)-^-\Tf-''-'^-HtdT 

■ J JlR2d 



N+d+l 
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for some constant Cq, > 0. This gives for j > J 
\ri''^{^^2Fj){r],x)\ 

(3.9) 



iv) 



-2(N+d+l) 



|(r/,x)p = (|r/|2+|xn"/^<2"/2|ry|"< (2rf)"/^ max \r 



\-f\<2(N+d+l) 

Now let V^ C M*^ be a closed conic set and let v? e C^(R'^'^) \ satisfy 

= r n (supp (f xv) = supp if X (unv), 

which by (13. 6p means that {ri,x) & V \0 ^ \x\ < \ri\. Let n E N\0 
and let {r],x) G V^ \ 0. Then we have 

'^ max 1""°"' 

a: \a\=n 

For j > J (13. 9p thus implies 
(3.10) 
sup \{r],x)\''\^{ip^2Fj){7],x)\<Cn sup snp\r]''^{ip^2Fj){r],x)\ 

(r;,x)gV {t],x)gV \a\=n 

<C„ sup sup a(r/)-2(^+'^+i)(a:)^+'^+^ 

{t],x)gV \a\=n 

<Cn sup SUpC„(77)-2(^+'^+l)(7;)^+'^+^ 
{ti,x)gV \a\=n 

for some C„ > 0, independently of j > J. This means that the second 
criterion (13. 4p ^ for the convergence (13. 7p has been proved. 

It remains to verify that ^2Fj,^2F G ^^(M^^). If (r/,x) ^ [/ and 
{r),x) 7^ then |x| < \ri\ and there is an open conic neighborhood U' 
containing {ri,x) of the form U' = {(^,x) : |x| < C*!?]!} for C > such 
that U n U' = ^. The estimate (I3.10p with V replaved by U' shows 
that WF{^2Fj) C r. Likewise WF{^2F) C L so we have proved 
^2Fj,^2F e^i^iR^'^). D 

Proposition 3.2. If F E C^o^iR'^'^) and t G R'^ then we have 

(3.11) ^{F{t, ■)) = {^2F){t, ■) zn ^'{R'^). 

Proof. The result (IXTTD says that ^ o i?^ = i?j o ^2 on C~ 



3 2d\ 



slowv 



where the equality is understood in ^'{R"^). Let U C M^'^\0 be defined 
by (13. 6p . Then F = R'^'^ xU is a closed set, conic in the second variable, 
and by (13. 2 p we have 

F n A^i = (t X R^) X f/ n (M'^ X 0) = 0. 

Thus the restriction (13. ip is a continuous map between distribution 
spaces according to (13. 5p . By Lemma [3?T] and (13. 5p we have {^2F)(t, ■) G 
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Let X e C^iW"), X > 0, x{x) = 1 for |x| < 1, Xj{x) = xi^/j) 
and Fj = Fxj ® Xj e C^{R'^'^). As before Fj ^ F in ^'(M^"!)^ 
and therefore also in ^'{R'^'^), as j — )■ oo. For any ip G C^(]R^'^) we 
have (p{Fj — F) = ip{xj ^ Xj ~ ^)F = ioi j sufficiently large, which 
together with Fj ^ F in ^'(M^^) j^piy that Fj -^ F in ^f (M^^)^ gj^ce 
Fj,F E ^f(M2'^) because iyF(F,) = WF{F) = 0. Thus (^ and 
(1^ implies that Fj{t, ■) ^ F(t, ■) in ^'{R'^). We have in fact 

(3.12) Fj{t,-)^F{t,-) in ^'(M"^). 

To see this, let ^ e ^(M^). Since |F(t,r)| < C(t)^(r)^ for some 
C, A^ > we have for any e > 

N I /^\N\ 



\{F,it,-),if{xk-m<2C{tr / {Ty'\cpiT)\dT<e uniformlyinj 

J|t|>A: 

for k > ki and /ci sufficiently large (t is fixed). Since F{t, ■) G ^'(M*^) 
and (pxk — !• v' in ^(M*^) as fc — )■ oo we have 

\{F{t,-),^{l-Xk))\<e iik>k, 

for A;2 sufficiently large. The latter two estimates now give for k > 
max(fci, ^2) 

|(F(t,-)-F,-(t,-),^)| 

< |(F(t, ■),^{1 - Xk))\ + |(i^(t, ■) - F,(t, ■),^Xfc)l + mit, ■),^{xk - 1))| 

< 36 

for j sufficiently large, since Fj{t,-) — )■ F{t,-) in i^'(]R'^). This proves 

Now we use Lemma [XT! The convergence (13.71) . (13. 5p and (13. 3p give 

{^2Fj)it, ■) — ^ i^2F)it, ■) in ^'(M'*) as j -^ 00. 

For ^ G C^(]R'^), this finally gives, using (13.121) and Fubini's theorem, 

(^(F(t,-)),^) = (F(t,-),^) 

= lim(F,(t,-),</^) 



lim / Fj{t,T)(p{T)dT 



lim / F,{t,T)( [ ^{Oe'^^-<dAdT 



j->oo 



lim / {^2Fj){t,0^{0d^ 



lim((^2F,)(t,-),^) 

J— >-oo 
((^2F)(t,-),^). 

D 
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If / e C,^^(M'^) then f ®J o k e CXo„(^^'^), so the definition of 
the Wigner distribution fll.ip combined with Proposition 13.21 gives the 
following byproduct. 

Corollary 3.3. If f e CX,^{R'^) and t G M"^ then we have 
^{f®JoK{t,-)) = Wf{t,-) zn^'(M'^). 
Finally Lemma 12.51 gives 

Corollary 3.4. If f e ^c?'(M"') and t G M^ then we have 
^{f®JoK{t,-)) = Wf{t,-) m^'{R'^). 

4. A WiGNER DISTRIBUTION MOMENT FORMULA FOR THE 
INSTANTANEOUS FREQUENCY 

Denote the modulus of2; = x + i?/GCbyr = \z\ and the argument 
(or phase) hj if = aigz. The polar-to-rectangular coordinate trans- 
formation on M^ ~ C is defined by {x,y) = g{r,ip) = {r cos ip,r sin (p). 
Denote the positive reals by IR+. Then each of the restrictions gi and 
g2, defined respectively by 

fi'l '■= fi'k+x(-7r,7r)5 5'2 '■= fi'k+x(0,27r)5 

is an analytic map with an analytic inverse, mapping surjectively on 
the open sets Ui and U2, respectively, defined by 

^1 : M+ X (-TT, tt) ^ C \ {(-00, 0] + iO} := f/i, 

^2 : M+ X (0, 27r) ^ C \ {[0, +cx)) + iO} := U2. 

The Jacobian of g is 



Dg 



which is invertible with inverse 




cos if —r sm ip 
sin ip r cos ip 



(4-1) ^j:— = -?— -^> ^5::^ = :;:j-:;:^' ^ + weUiUU2. 



/^ N_i _ 1 / rcos(y9 rsinv9 \ _ 1 / x^/x'^ + y"^ y\/x'^ + y"^ 
r\-sm.ip cosip ) x"^ + y'^ \ -y x 

provided r > 0. The inverse function theorem applied to each of the 
restrictions gi and g2 thus gives 

daigz —y daxgz x 

dx x"^ + y"^^ dy x^ + y"^ 

(Note that gi^{z) and g2^{z) have arguments that differ by 2ii for 
Im z < 0, but this does not affect the partial derivatives of arg^;.) 

Let fit) = u{t) + iv{t) be a function / : M"^ 1— )■ C. If / is continuous 
then [/j,i := {t e W^ : f{t) G Ui) and Uf^2 ■= {t e R'^ : f{t) G 
U2} are open sets, and it follows from above that each of arg/(t) = 
igi' o /(t))2 : f//,i ^ M and arg/(t) = (g^' o f{t))2 : t//,2 ^ M 
are continuous. Here (g^^ o /(t))2 means the second component of 
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gi\fit)). life C^R"^) then, likewise, each of arg /(t) : f//,i ^ R and 
arg/(t) : Uf^2 ^-^ R sue differentiable. From (14. ip and the chain rule we 
obtain the partial derivative with respect to tj, in both cases, as 



(4.2) 



daxg fit) daxgz du{t) daxgz dv{t) 

=f{t) dtj 



dtj dx 



=/(t) dtj dy 

u{t)dtAt) - v(t)dt,u(t) 



u\t)+v^{t) 
for t e f//,i U f//,2 = {teR'^: fit) ^ 0} and I < j < d. Thus 

^--^/")- '""!;i:!;::i;r' ^<'>^°- 

The instantaneous frequency [IJ of f{t) is defined as the normalized 
gradient (27r)~^Varg/(t) with domain {t : f{t) ^ 0}. For a character 
^2m£,-t ^j^]-^ frequency ^ the instantaneous frequency is thus ^ constantly, 
which means that the instantaneous frequency is a generalization of the 
concept constant (global) frequency. 

There is a connection between the instantaneous frequency of a suffi- 
ciently smooth and decaying function and its Wigner distribution [Hll] . 
In fact, the instantaneous frequency can be written, for fixed t e R'^, 
as a normalized frequency moment of order one of the Wigner distri- 
bution. (In the engineering literature a heuristic version of this result 
is well known [T].) 

Proposition 4.1. Suppose e > and f G H2+^+'^(R'^). Then for any 
t eR'^ such that f{t) 7^ we have 

(4.3) fva../W^ ^'W^ 

Proof. First we note that / G if 2+i+<^ C L^ and (12. 3 p imply that 
Wf G Co(M^'^). Hence the restriction of the Wigner distribution Wf H- 
Wf{t, ■) is a well defined continuous map C^iR?'^) t-)- Co(M'^) for any 
t gM'^. 

Let t G M*^ be fixed arbitrary such that f{t) 7^ 0. The assumption 
/ G iif 2+^+'^(M'^) and the Sobolev embedding theorem imply that 

(4.4) /, 9,/ G Co(M"), l<3<d. 
Cauchy-Schwarz and f E L"^ give 



(4.5) g:=f{t + -/2)f{t--/2)eL\R''). 

The assumption / G if 2^"'^^^(]R°') also gives 
(4.6) 

/ mmdi <( [ {o'^-^^'^'^\m\'d^) ( [ (o-'-'^dA <oo. 

jRd. XJRd. J XJRd J 
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Hence / G L\R'^) and since ^(/ (t + ■/2)) (0 = 2'^Mtf{2^) w e have 
^(/(t + 72)) G Li(M'^). Thus ^(/(t--/2)) = ^(/(t + 72)) imphes 
that 

(4.7) ^(? = ^(/(t + 72)) * ^(/(t + 72)) e L'iR"). 

By dOD, 033]) and (iZD we have 5- G (ConLin^Li)(]R'*), which means 
that the conditions for Fourier's inversion formula to hold pointwise [H] 
for the function g are satisfied. In particular 

(4.8) \f{t)\' = g{0)= [ g{Od^= [ Wj{t,Od^, 

in the last step using the definition (11.11) . 

Concerning the numerator in (14.31) we obtain using integration by 
parts and the fact (14.41) that / vanishes at infinity 

(4.9) ^r r 

■%i^f{t + r/2)f{t-r/2))dT 



27ri ' 



for I < j < d. Let us study the function gj := djf{t + ■/2)f{t — ■/2) 
and 

hAr) ■■= dr, {fit + r/2)f{t-r/2)) = ^{g,{r) - ^;{=t)) 

ioTl<3 <d fixed. We have ^^/(O ^ /.^(M"') and thus djf G L^R"^). 
Since / G L^(M'^) the Cauchy-Schwarz inequality gives gj G L^{R'^). 
From above we know that J^{f{t — ■/2)) G L^(]R'^) and likewise we 
have 

^{d,f{t + -/2)m = 2V-^«-*5/(20 = 2''2me'^'^-%,m) e L\R') 
because of (14. 6p . Thus 



g, = ^{d,f{t + 72)) * ^{f{t - 72)) G L'iR''). 

Invoking flO|) we have proved that gj G (Cq fl L^ fl ^L^){R'^) which 
means that Fourier's inversion formula holds for gj and thus hj. Inte- 
gration of (14. 9 p gives, denoting f(t) = u(t) + iv(t), 

1 f ^uu^.c-hM 



[ ^^Wf{t,Od^=^ I ^h^Od^ 



2ni 



9M-9M 



(4.10) 4:TTi 



= J- [d^mm - d, fit) fit) 

= ^ iuit)d,vit) - d,uit)vit)) . 
Finally (gSD follows from a combination of (giS]), (iSD and KTO\\ . D 
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Next we will prove a version of Proposition 14.11 for functions / G 
^(S"{M.'^). Note that this assumption is neither a generalization nor a 
special case of the assumptions of Proposition I4.1[ 

Proposition 4.2. If f e ^^'(M"^) then 

(4.U) _Lvarg/W^<^|^. /(t) # 0. 

Proof. Let t G M'^ be fixed such that /(t) ^ 0, and let f{t) = u{t)+iv{t). 

Set 

g,:=f(t + -/2)f{t--/2). 

By Corollary [331 we have ^gt = Wf(t,-) in ^'(R'^). The assump- 
tion / G t^<f"(]R'^) together with the Paley- Wiener-Schwartz theorem 
implies that / extends to an entire function on C^ and there exist 
N,A,C >0 such that 

\f{x + iy)\<C{x)^e^\y\ x,yeR'^. 

Hence the function gt extends to an entire function on C^ with bound 



\gt{z)\ = \fit + z/2)fit-z/2)\ 

< C{t + Re z/2f{t - Re z/2fe^^^^'^ "/^l 
<C(t)2^(Rez)2^e^|i-^l. 

By the Paley-Wiener-Schwartz theorem ^gt G S"{M.'^) and ^gt is 
supported in a fixed compact set independent of t G M''. Since ^gt = 
Wflt, ■) in S>'(M.'^) we have Wf{t, ■) G ^(M'^). 

Because Wf(t,-) has compact support, {^gt,l) and {^gt,i) are 
well defined, where {^gt^i) means the vector {{^gt,^j))'j=i and ^j : 
M*^ I— >■ M is the jth coordinate function. The inverse Fourier transform 
of a distribution u G S"(M.''-) is the entire function x H- {u, e^'^*^') [8], 
Theorem 7.1.14] so we have 

gtiT) = ^-'^gtiT) = {^g,,e'-n 
and in particular 

(4.12) {Wfit, ■), 1) = {^gt, 1) = ^*(0) = |/(t)P = u'it) + v\t). 

In a similar way we have for 1 < j < d 

1 

~ 27n 
(4-13) = ^^ 



-(^(%,),1) 



^ (a,/(t)/(t) - /(t)9,/(t) 



1 

2^ 



— (^W^'.^W - dAt)<t)) 
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Now IKWf follows from IKWf . (KT^ and (i2D. □ 



Finally we give a generalization of Proposition 14.21 
Proposition 4.3. Let 

f{t) = exp{7rit ■ At)h{t) 
where A G W^^'^ is symmetric, and h G ^<S"{M.'^). Then 

(4,4) -Lv..g/W.iS|M. m^o. 

Proof. Let t G M*^ be fixed and satisfy f{t) ^ 0. We have 
arg /(t) = arg h{t) + ixt ■ At, 

(4.15) V arg f{t) = V arg h{t) + 2nAt, 

and 

Wfit,0 = W,it,^-At). 
This yields 

(4.16) {Wfit,-)A) = {W,it,-),l) 

and 

{Wf{t,.),Q = {W,{t,--At),Q 
^ ■ ' ={W,{t,-),Q + {At),{W,{t,.),l). 

Finally dUSD, ( KTf} . Proposition S2] and flJAHj) give 

{Wfit,.),Q (iy,(t,-),o) + (At),(iy,(t,-),i) 
(iy/(t,-),i) (w^/.(t,-),i) 



{Wh{t,-),1) 



+ Wj 



— Varg/i(t) + At 
2tt ,^ 

7^Varg/(t) 
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